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Abstract. A renormalization scheme for interacting fermionic systems is presented where the renormaliza- 
tion is carried out in terms of the fermionic degrees of freedom. The scheme is based on continuous unitary 
transformations of the hamiltonian which stays hermitian throughout the renormalization flow, whereby 
any frequency dependence is avoided. The approach is illustrated in detail for a model of spinless fermions 
with nearest neighbour repulsion in one dimension. Even though the fermionic degrees of freedom do not 
provide an easy starting point in one dimension favorable results are obtained which agree well with the 
exact findings based on Bethe ansatz. 
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1 Introduction 

With the discovery of high temperature superconductiv- 
ity in layered cuprates of perovskite type 1 and the sub- 
sequent most intensive theoretical considerations of this 
intriguing phenomenon (see e.g. Refs. 0E]) it has become 
apparent that a reliable approach to strongly interacting 
fermionic systems is lacking. So the efforts to formulate 
the very successful concept of renormalization 0] also for 
extended interacting fermionic systems have been intensi- 
fied considerably during the last decade, see e.g. Ref. [S]. 
By now the literature on this approach is so wide that it is 
not possible to provide an exhaustive list. This underlines 
the importance that is accorded to this topic. 

So far a number of renormalizing schemes has been 
applied to two-dimensional Hubbard-type models which 
are relevant for high temperature superconductivity 
EE]- These schemes rely conceptually on diagrammatic 
perturbation theory in the interaction strength. They are 
non-perturbative in the sense that infinite orders of the 
interaction are kept; the necessary truncation concerns 
terms of a certain structure, e.g. six points correlations, 
not terms of a certain order in the interaction. The Fermi 
surface is discretized and the various scattering couplings 
across the Fermi sea are suitably parametrized. It is pos- 
sible to detect whether or not the renormalized couplings 
decrease of diverge in the course of the flow. In this way, 
robust evidence for the occurrence of <i-wave supercon- 
ductivity was found. The corresponding couplings diverge 
for certain values of doping and interaction. So far, how- 
ever, no calculations exist for dynamical correlations like 
the ARPES response at finite energies and wave vectors 
(measured relative to the Fermi surface) as required for 



the understanding of the experimental findings. This is 
due to the great complexity of the problem which requires 
- among other difficulties - to follow the flow of the ob- 
servables as well, see for instance the appendix in Ref. [B] • 
For an impurity in a spinless Luttinger liquid a spectral 
density at all energies was computed by a one-particle ir- 
reducible renormalization approach neglecting, however, 
the renormalization of the two-particle vertex [5]. 

Starting from a flow equation approach as proposed 
by Wegner JHj and quite similarly by Glazeck and Wil- 
son |11I12| a renormalizing scheme based on continuous 
unitary transformations (CUTs) has been proposed |13| . 
By an appropriately chosen unitary transformation an ef- 
fective hamiltonian H e g is obtained. The transformation 
is tuned in such a way that H c q conserves the number of 
quasi-particles. Due to this property the calculation of dy- 
namical correlation function for all energies becomes pos- 
sible as was shown for spin ladders |14j . For this reason, we 
consider the renormalizing CUT approach to have a par- 
ticularly great potential. This expectation is supported 
decisively by recent work on quantum chemical systems 
[*Po| where White could show that the numerical applica- 
tion of a continuous unitary transformation similar to the 
one used in Ref. ^Hj an d nere leads to excellent results. 

It is the aim of the present work to explain the techni- 
cal details of the calculations announced in Ref. To 
this end, intermediate results will be shown. We hope that 
the available data will make it possible to conceive also 
analytical treatments which capture the essential physics. 
The medium-term objective is to generalize the renor- 
malizing CUT approach to more realistic models. A de- 
manding challenge is to treat the dimensional crossover 
between one- and two-dimensional interacting fermionic 
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models. This is of great experimental relevance since the 
physical systems existing in nature are at best quasi-one- 
dimcnsional, i.e. strongly anisotropic, so that the higher 
dimensionality enters always at a certain stage. 

Two dimensions are the most demanding case for the 
theoretical description of strong correlations. In three di- 
mensions the powerful Fermi liquid theory is well estab- 
lished which is based on the observation that the quasi- 
particles as such yield a good description of the low-lying 
excitations (see e.g. The interaction between the quasi- 
particles is not essential. In one dimension on the other 
hand, the collective plasmon modes dominate the low- 
energy physics completely so that the fermionic hamil- 
tonian can be mapped to a bosonic one representing so- 
called Luttinger liquids (see e.g. [T7IT%] ). This phenomenon 
can be seen as a binding (or anti-binding) of a pair con- 
sisting of a hole and a fermion. This implies that it is a 
signature of a dominating interaction between the quasi- 
particles. This fact is commonly interpreted as the failure 
of a description of the energetically low-lying physics in 
terms of quasi-particlcs. 

Considering two dimensional strongly interacting sys- 
tems it is shown that they are generically Fermi liquids 
in the weak coupling regime |19j . So, qualitatively, two di- 
mensional systems are similar to three dimensional ones in 
the weak coupling limit. But at strong coupling this is no 
longer true. For given generic hopping element t and in- 
teraction strength U the lower dimensional system is more 
influenced by strong correlations than the corresponding 
higher dimensional one. This stems from two effects: (i) 
The band width being proportional to the coordination 
number is higher in the higher dimensional case, (ii) If 
collective modes, damped or undamped, are formed their 
density of states (DOS) at low energies is higher in lower 
dimensions. For instance, the DOS p(u>) of linear dispers- 
ing modes u> oc |k| behaves like p(ui) cx uj d ^ 1 with dimen- 
sion d. 

In the above sense, two dimensional strongly interact- 
ing systems represent an intermediate situation. Generi- 
cally, neither the interactions between the quasi-particlcs 
can be sufficiently described by a Landau function as in 
three dimensions, nor is the physics completely dominated 
by collective modes formed from bound particle-hole pairs. 
Hence, in two dimensions one has to have a theoretical 
tool which is able to reconcile both main features: col- 
lective modes occur and they are important, but they do 
not exhaust all degrees of freedom. There are also quasi- 
particles. But their interaction is very important. 

The above considerations are the motivation to show 
here that it is possible to use continuous unitary trans- 
formations and quasi-particle description for one dimen- 
sional systems to recover the known results. In particular, 
we will look at the momentum distribution in the ground 
state which differs significantly between Fermi liquids and 
Luttinger liquids. In Fermi liquids a jump by Zk F , the 
quasi-particle weight, occurs at the Fermi level whereas in 
Luttinger liquids only a power-law behaviour occurs |17l 
I18| . Our investigation is intended to be a test case for 
the method, not as a means to obtain new and so far un- 



known data. Evidence is provided that in one dimension 
the physics has not turned bosonic out of the blue but that 
a description in terms of fermionic quasi-particles is still 
reasonable even though a description in terms of bosons 
is easier. In view of the medium-term aim to describe the 
dimensional crossover we consider the fermionic approach 
to be a necessary prerequisite. 

The paper is set-up as follows. After this Introduc- 
tion the method is described in Sect. 2. Also the model to 
which the continuous unitary transformation is applied is 
given in detail. In Sect. 3 the numerical results are shown 
and discussed. The comprehensive Discussion concludes 
the article in Sect. 4. 



2 Method and Model 

In general, we consider a translationally invariant system 
of N interacting fermions 



H = NE + £ «: C ' C , 



kqp : C k+q C k-q C k-p C k+p ' ' W 



kqp 



Here c' k (c.) creates (annihilates) a fermion at wave vec- 
tor k in momentum space. Note that the parametrization 
of the scattering processes is not the conventional one. 
Our choice, however, is more apt to represent the inher- 
ent symmetries (see below). The fermions appearing are 
considered to be spinless, i.e. there is at maximum one 
per site and no spin index occurs. The colons : . . . : de- 
note normal ordering with respect to the non-interacting 
Fermi sea. This commonly used classification makes it eas- 
ier to trace the effect of the individual terms (cf. appendix 
A). The function r^qp is the vertex function embodying 
the amplitudes for all possible scattering processes. 



2.1 Method 

The method which we use to analyse the hamiltonian 
is the method of continuous unitary transformations 
(CUTs) based on the flow equation approach proposed by 
Wegner . The unitary transformation U (£) is parametrized 
by the continuous variable £ ranging between and oo. A 
given hamiltonian H(0) is continuously mapped onto an 
effective model H c ff as £ is taken from zero to infinity. The 
continuous unitary transformation is defined locally in £ 
by the antihermitian generator rj(£) 



-H(£) = [r ] (£),H(£)} 



(2) 



Of course, all other observables O have to be subject to 
the same transformations 



dl 



o{t) = W),o{i)] 



(3) 



CP. Heidbrink, G.S. Uhrig: Renormalization by Continuous Unitary Transformations 



3 




Fig. 1. Hamiltonian with a block band structure transformed 
into an effective block diagonal hamiltonian. Each black block 
stands for the action of the hamiltonian in a subspace of the 
Hilbert space with a given number of excitations. Let us as- 
sume that the uppermost block stands for the vacuum (no 
excitation); the next lower block acts on states with one sin- 
gle excitation; the next lower block acts on states with two 
excitations and so on. In the example depicted the number of 
excitations can change at most by two on a single application 
of the hamiltonian. This is the block band structure. 

since the expectation values and correlations shall not be 
altered by the transformation. 

The main task is to determine rj(£) — i](H(£)) in a 
way that brings the hamiltonian systematically closer to 
a simpler structure. Once the generator is chosen and the 
commutator calculated, one obtains a high dimensional 
set of coupled ordinary differential equations. These have 
to be solved in the limit £ — > oo. 

For the choice of r](£) we focus on the case where the 
original hamiltonian H has a so called block band struc- 
ture with respect to a certain counting operator Q, i.e. an 
operator with a spectrum of non-negative integers. The 
operator Q shall be the operator counting the number of 
excitations present in the system. The ground state |0) of 
Q, the so-called quasi-particle vacuum, is the state with- 
out any excitation so that 

Q\0) = (4) 

holds. The block band structure occurs if the whole hamil- 
tonian H changes the counting operator Q at maximum 
by a certain value M < oo. This means that H links two 
eigen states \i) and \j) of Q only if their eigen values qi 
and qj differ at most by M: |g, — qj\ < M. In Fig. Q the 
case M = 2 is illustrated. 

In view of systems of interacting fermions we choose 
Q to be the operator counting the number of quasi-particles, 
i.e. the number of particles above the Fermi level plus the 
number of holes below 

Q = 5>gn(|fc| - fc F ) : cK : ■ (5) 
k 

This choice implies that we intend to use ordinary quasi- 
particles as the elementary excitations. Generally, we like 
to stress that the choice of Q requires some physical intu- 
ition or a certain presumption about the problem under 
study. The choice of Q can be considered as the choice of 
a starting point. Depending on the quality of this choice 
subsequent approximations will be more or less reliable. 
This point will be discussed below in more detail. 



First, however, we turn to the choice of the generator. 
For band-diagonal matrices Mielke proposed an ansatz 
|20j that was generalized to many-body problems with 
block-band structures by Knetter and Uhrig (23- The an- 
tihermitian generator 7](i) is chosen such that its matrix 
elements are given in an eigen basis of Q by 

r Hj (t) = S ign(q i (£)-q j (l))H ij (£) . (6) 

For the interacting fermions as in the conventional 
occupation number basis is appropriate. 

Since in general the eigen space for a given number of 
excitations qi has a finite dimension we use the convention 
that Aij is not only a single matrix element but the whole 
submatrix of A which connects the eigen space belonging 
to qj (domain) to the eigen space belonging to qi (co- 
domain) . With this convention Eq. [H| becomes a matrix 
equation. Inserting Eq. into the general flow equation 
yields 

-TjHij = -signfe - qj)(HuHij - //,,//, , ) 

+ ^2 ( si S n (* ~ Ik) + sign(<?j - q k )) H lk H kj (7) 

k^ij 

which is also a matrix equation, i.e. the sequence of the 
matrices in the products matters. One can show [2THl2*T] 
that 

— preserves the block band structure and 

— leads to a block diagonal H c fi where qi =/= qj =>■ H c g ^ = 
0. The condition necessary for these conclusions is that 
the spectrum of the hamiltonian is bounded from be- 
low which constitutes a natural assumption for physi- 
cal systems. 

The block diagonality of the effective model is equiva- 
lent to the commutation of the effective model with Q 
[Q,H eS ]=Q. 

Furthermore, we show that the ground state of the 
system is mapped onto the state with no elementary ex- 
citations, i.e. the vacuum |0) without any excited quasi- 
particles, in the course of the transformation. Let us as- 
sume that the ground state is not degenerate. Without 
loss of generality the indexing is chosen such that i = 
refers to the vacuum with qo = 0. For j > = we obtain 
from 

f l H i0 = -(H ii H i0 -H i0 E) 

+ ( si g n fe ~ 9k) - 1)) H lk H k0 , (8) 

k=£iO 

where E := (0|i?|0) = Hqq is the energy of the vacuum. 
H i0 is a vector that connects the vacuum with other states. 
For large £ the flow has led the system already close to 
block diagonality |20U21j . Then the matrices Hi k linking 
blocks of different number of elementary excitations are 
small quantities so that Eq. (JSJ is dominated by the first 
term on the right hand side. The products Hi k H k Q in the 
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sum are smaller by one order in the off-diagonal matrices. 
So the asymptotic behaviour is given by 

^H lQ ~ (E-H u )H i0 . (9) 

Since Hiq has to tend to zero due to the general conver- 
gence |2T)ll2"l"] we deduce from that E—Ha < for large 
£. This implies that E is lower in energy than any state dif- 
ferent from the vacuum which is linked to the vacuum by a 
non-zero element HiQ. We interprete the physical content 
of this result in the following way. In a block band diag- 
onal hamiltonian as illustrated on the left side in Fig. ^ 
the vacuum is linked to states with a certain finite num- 
ber (bounded by M) of elementary excitations. From 10 
we learn that the vacuum is lower in energy than those 
states. Generically, this implies that also all other states 
which contain an unrestricted number of elementary exci- 
tations lie higher in energy. We conclude that the vacuum 
is indeed the ground state unless a phase transition takes 
place. 

To illustrate the latter statement in more detail let us 
think about a hamiltonian which is controlled by an inter- 
action parameter U such that the vacuum is the ground 
state for U = without any transformation. Then the 
situation changes gradually and smoothly as U is turned 
on so that the ground state is mapped onto the vacuum. 
This is true till a singularity occurs, that is a phase tran- 
sition. A second order phase transition would be signaled 
by the softening - vanishing of the eigen energy - of one 
of the excited states (for an example see the single triplet 
excitation in Ref. [221 )• Technically, this implies that the 
convergence of HiQ as given by Eq. © breaks down since 
the expression in the bracket on the right hand side van- 
ishes. This reflects the physical fact that our approach will 
not work beyond the phase transition since there the orig- 
inal vacuum is no longer a good reference state for the 
true ground state. But it is possible to use the approach 
until the phase transition is reached and the second order 
phase transition can be detected by the vanishing of the 
energy of one of the excited states. 

First order transitions spoil also the applicability of a 
smooth mapping. But they cannot be detected locally. A 
first order jump occurs when a completely different state 
comes down in energy. Such a completely different state is 
built from a macroscopic number of elementary excitations 
and will in general not be connected to the vacuum by a 
finite element HiQ. This means that the smooth mapping 
constructed by the CUT may still work even though a first 
order transition has occurred simply because there is no 
local connection of the local energy minimum, which is 
mapped to the vacuum, to the global one. This concludes 
the general considerations about the mapping generated 
by the CUT between the ground state and the vacuum of 
elementary excitations. 

Now we return to the model of interacting fcrmions in 
|QJ. Inspection shows that this model has a block band 
structure with respect to the counting operator (J5J of 
quasi-particles. If a particle from below the Fermi level 
is scattered to a state above the Fermi level two elemen- 
tary excitations are created: one hole and one particle. 



The inverse process corresponds to the decrement of the 
number of quasi-particles by two. If the particle is scat- 
tered from below the Fermi level to below the Fermi level 
the number of quasi-particles remains constant. The same 
is true if the particle is taken from above the Fermi level 
to another state above this level. If we consider at maxi- 
mum scattering terms built from four fermionic operators, 
as it is done in Eq. Q , then two fermions are scattered so 
that the possible changes in the number of quasi-particles 
are 0,±2,±4. So there is a natural upper bound M = 4 
and the system displays indeed a block band structure. 
Note that all the statements so far were independent of 
the dimension of the model. 



2.2 Model 

The explicit model we consider here is a system of one- 
dimensional spinless fermions. Haldane used it to explain 
the concept of Luttinger liquids [231 ■ Physical realizations 
one may think of are either completely polarized electrons 
or anisotropic spin chains which can be rigorously mapped 
onto spinless fermions by means of the Jordan- Wigner 
transformation |24j . A completely different application is 
the description of vicinal surfaces where the hard-core re- 
pulsion between different steps is taken into account by 
passing to fermions |25| . 

The model is a tight-binding model where the particles 
are distributed over a one-dimensional chain of N sites 
with antiperiodic boundary conditions at half filling, i.e. 
on average each site is occupied by half a fermion. The 
fermions can hop to nearest neighbour sites and there is 
a repulsive interaction V between two adjacent fermions 

N 

H = 12 [j(4fi«* + h - c -) + V " I) ("< " I) . 

i=l 

(10) 

where a\ {a i ) creates (annihilates) a fermion on site i in 
real space. This model is exactly solvable in the form of 
an anisotropic spin model. The solution is due to Bcthc 
and to Yang and Yang |27ll2^| . More results were 
found later 29 ,30 ,31, 32 23]. In the present context it is 
important to know that the system is metallic for not too 
large couplings, namely for V < 1. In this region it is 
the simplest realistic tight-binding model displaying Lut- 
tinger liquid behaviour which is why we use it as our test 
case. At the value V = 1 a continuous phase transition 
takes place into a charge density wave where the fermion 
density is staggered. Thus the translation symmetry of 
the system is spontaneously broken. This phase displays 
also a charge gap and represents hence an insulator. In 
the present work, however, our interest is focused on the 
metallic phase. 

Representation Let us first clarify the notation that we 
use in the following. Transforming the real-space represen- 
tation l|10|) into momentum-space and normal-ordering of 
the hamiltonian with respect to the non-interacting Fermi 
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sea yields a hamiltonian of the form in Eq. with the 
coefficients 



E 

£k 

■^kqp 



-(l + ^)cosfc 
V sin q sinp . 



(11) 
(12) 
(13) 



in the thermodynamic limit N — > oo. For finite system 
sizes TV < oo as used below the dispersion reads 



£fc 



2V_ 
~~N 



1 H cos k' cos k 



E 

fe'i<fc 



(14) 



The terms - and s- in J5 correspond to the ground state 
energy of the non-interacting system and to the Fock term, 
respectively. We are not interested in these leading order 
effects but will concentrate on the correlation effects as 
they appear in the subsequent orders. 

Based on the parametrization chosen in Eq. Jy the 
symmetries and notation properties of the system manifest 
themselves in a very concise way in the vertex function 

^kqp 

1. hermitecity leads to 



Fkqp — ^kpq 7 

2. inversion symmetry leads to 

-^kqp — F—k — q- 

3. particle- hole symmetry leads to 



Fkqp 



k+7rqp 



= r, 



(15) 



(16) 



(17) 



A swap of two neighboured fcrmionic operators in a normal- 
ordered product leads only to an additional minus sign. 
Thus there is a redundancy in the notation : c\Ci ...:= — : 
C2C1 . . . : which implies that Fk q p is not uniquely defined 
by Eq. . This caveat can be remedied by the additional 
requirement of the 

4. notation symmetry 



r^qp = — -T/ 



k-qp 



= -fib 



kq-p 



(18) 



The notation symmetry implies that the momentum de- 
pendence of the vertex function is the one given in (|13|) . 
Note that in our notation the vertex function r^qp van- 
ishes automatically if fermions are scattered from the Fermi 
points to the Fermi points (k = ±ir/2 and q,p £ {0, ±7r}). 
This important fact leads to the property that the inter- 
action is not a relevant perturbation but only a marginal 
one. The CDW does not occur at arbitrarily small inter- 
action but only beyond a finite threshold. If the scattering 
is denoted in the usual way a much more elaborate rea- 
soning computing the contribution of the Cooper pair and 
the zero sound channel is required to yield the same result 

ma. 

Exploiting the above symmetries the number of non- 
zero scattering amplitudes that has to be dealt with is 
reduced from N 3 to approximately 



Exact Results Here we report some exact results to which 
we will compare our findings. 

1. The ground state energy per site in the metallic phase 
V < 1 reads [23 



E 
N 



(V) 



COS fl 



sin 2 /i 



dx 



2cosh7ra:(cosh2x/i— cos/x) ' 

— oo 

(19) 

where V — cos/i. Ea. l|19|l is dominated by its linear 
part — (— + j-) (cf. Eq. lfTT|) '). which we will substract 
to focus on the correlation part of the energy. 

2. We deduce from the dispersion of the anisotropic Heisen- 
berg model |29j the dispersion of the effective fermionic 
model 

£k = - f^-jcosk , (20) 

which we expect after the CUT. In particular, we use 
the value of the Fermi velocity as done previously 
|32j . This quantity is also dominated by a rather trivial 
linear term 1 + 2Z (cf. Eqs. (|12I14I> ) which we will 
subtract in order to focus on the correlation effects. 

3. The momentum distribution n(k) in the ground state 
cannot be computed directly with Bethe ansatz. But 
it is possible to deduce the asymptotic behaviour for 
momenta close to the Fermi wave vector fcp based on 
the representation of the model in terms of bosonic 
degrees of freedom 32 . In the momentum distribution 
a characteristic signature of Luttinger liquid behaviour 
shows up. No finite jump in n(k) exists at k = ±fcp but 
a power law singularity with non-universal exponent 
appears. This singularity is of the form (for more 
details, see e.g. |3*5HT7] ') 



n{k ~ k-p) 



Cisign(fc — kp)\k — ky\ 



with 



a(V) 



T)(V) 



- 1 



(21) 

(22) 
(23) 



7T — arccos V 
and an undetermined constant C\. 

2.3 The Continuous Unitary Transformation 

In the subsection on the method \l. II it was explained that 
the block band structure of the system is preserved under 
the chosen CUT. This means that the number of quasi- 
particles is altered at maximum by 4. It must be empha- 
sized, however, that this does not imply that only terms 
made from four fermionic operators, so-called 2-particle 
operators 1 , occur in the hamiltonian. Even though the ini- 
tial hamiltonian contains only 0-particle (constants), 

1 The name stems from the property that the effect of this 
operator is to re-distribute two real (not quasi) particles. This 
is due to the conservation of the number of real particles. 
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1-particle (of the form : c'c :) and 2-particle terms (of the 
form : c'c'cc :) the transformation may and will generate 
also 3- (and more) particle terms. But these terms have 
to fulfill the condition that they change the number of 
quasi-particles by no more than ±4. 

Let (£) be the normal-ordered j-particle term in 
H(£) which changes the number of quasi-particles by i. So 
the general structure during the transformation is 



H{£) =T (0) (£)+T Q (1) (£) + 



T 



0) 



+ 2 



T, 



0) 



T 



T 



U) 



J>2 



and 



]rsign(z)7; ( %) 

i>2 



{£) (24) 



(25) 



Note that the 1-particle term cannot change the number of 
quasi-particles since the conservation of momentum does 
not allow to shift the particle in momentum-space. 

Some general analysis is possible. One can study which 
combinations of Tr ' and Tr, -* in [rj(£), H(£)] influence 

ffTy' . First, the chan ges in the number of quasi-particles 
is additive i = i' + i" • Second, the maximum value of j is 



Jn 



1 



(26) 



due to the properties of the commutator between products 
of fermionic operators. Third, the normal ordering yields 
also a lower bound for j because the number of possible 
contractions is restricted. Contractions in a product made 
from normal-ordered factors are possible only between the 
fermionic operators from different normal-ordered factors, 
see Appendix A. So we are led to [3*6*| 



Jmin — 

1 j(f+i') + (j"-i") + \(f-i 

2 I (f - i') + (j" + i") + \(f + V) - {]" - i") 



(27) 



(j" + »") 



It is also possible to derive a set of abstract differ- 
ential equations for the coefficients of all possible terms. 
But its structure is quite complicated j36| . So a solution of 
the complete CUTs seems to be hardly accessible and we 
refrain from persuing this route further and turn to a nu- 
merical treatment of finite-size systems with N sites. Still 
further approximations are necessary because the number 
of different j-particle processes grows generically as 2V J-— 1 . 



2.4 Self-Similar or Renormalizing Approximation 

We analyse the flow of all terms of the form 

T^ ] {£) =E{i) (28) 
4 l) {t)= £ e k (£):clc k : (29) 

fce[0,27r) 

t} 2 \£)= W) ■■ A + A-, c k-v c k +P ■■ (3°) 

fce[o,7r) 

gpE[0,27r) 



with i G {0, ±2, ±4}. Terms involving interactions deal- 
ing with more than two particles T-^ >2 ^ are neglected. 
So we proceed as follows. The sum (|24f) of the terms in 
Eqs. (|28I29I30|I represents the hamiltonian at a given value 
of £. The corresponding generator r\ in Eq. I)25|l contains 
the same terms. This ansatz is inserted at the right hand 
side of the general flow equation @. The commutator is 
computed and the result is normal-ordered (cf. appendix 
A). The 3-particle terms are omitted and the coefficients of 
the 0-, 1- and 2-particle terms are compared. In this way, a 
high dimensional set of ordinary differential equations for 
the coefficients of the terms in Eqs. I|28I29I30[1 is obtained. 
These equations are given in appendix B. At £ = oo the 
parts of the third term l|3U|l which alter the number of 
quasi-particles will have vanished so that only the i = 
part remains (for illustration see Fig. EJ) . At this stage, i.e. 

at the end of the transformation, Tg°' ) is the ground state 

energy, the I-particle dispersion and represents 
the interaction of two quasi-particles. The latter does not 
need to be small. 

Since the given structure of the initial hamiltonian 
H(£ — 0) is preserved we call this approximation "self- 
similar" in the spirit of the work by Glazeck and Wilson 
The naming "renormalizing" is based on three facts. 
First, on the technical level the coefficients appearing in 
the initial hamiltonian are changed, i.e. renormalized, in 
the course of the transformation. Second, the procedure 
is non-perturbative since terms are omitted not because 
they are of a certain order in the initial interaction V 
but because of their structure being 3- (or more) particle 
terms. This implies that the couplings kept acquire infi- 
nite orders in V. Third, the generator (|6I25|I which we use 
here leads to a smooth exponential cutoff exp(— |Z\_E|^) 
of the matrix elements connecting states of different en- 
ergy (energy difference AE) |2UII21) . Thus matrix elements 
between energetically distant states are suppressed much 
more rapidly than those which are energetically very close 
to each other. This is similar to what is done in Wilson's 
renormalization 0] where the degrees of freedom at large 
energies are integrated out first. 

We illustrate the exponential cutoff for matrix ele- 
ments connecting to the ground state. Their asymptotic 
behaviour on £ — > oo is governed by Eq. ((5J)- Let us as- 
sume that the non-diagonal elements HiQ are already very 
small. Then the diagonal energies like E and Ha devi- 
ate from their asymptotic values only quadratically in 
the non-diagonal elements as results from second order 
perturbation theory. So in leading order the diagonal ele- 
ments can be considered constant in £. Then Eq. @ yields 
HiQ oc exp(— (Ha — E)£) as stated before. 

How can the restriction to the terms in Eqs. 



be justified? One argument results from considering an ex- 
pansion in the interaction V . The 3-particle interactions 
neglected are generated by a commutator of two normal- 
ordered 2-particle terms which are of the order of V so that 
they are of the order V 2 . For our purposes it is important 
to know in which order the terms kept are influenced by 
the neglect of the 3-particle terms. The 3-particle terms 
can have an influence on the terms kept only if they are 
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commuted again with at least a 2-particle term of the or- 
der V or higher. Hence the neglect of the 3-particle terms 
introduces deviations only in order V 3 . This holds for the 
1- and 2-particle terms (I29I30[1 . The O-particle term, which 
becomes at I = oo the ground state energy, is not changed 
by the commutator of a 2-particle and a 3-particle term 
since the generated terms cannot be contracted completely 
as stated by Eq. I|27|) . So the deviation in the O-particle 
term engendered by the neglect of the 3-particle terms 
is at worst of order V . Thus the approximation can be 
justified for low values of V. 

Another argument, which is more general than a power 
counting in the interaction, comes from the structure of 
the terms neglected. Let us focus on the 3-particle terms. 
Since they may change the number of quasi-particles at 
most by 4 they contain at least one annihilator of a quasi- 
particle. Since they are normal-ordered this annihilator 
appears rightmost. Hence such a term is active only if ap- 
plied to a state which contains already some excitations. 
Hence the approximation chosen is justified if the system 
can be considered a dilute gas of quasi-particles irrespec- 
tive of the strength of the interaction between the quasi- 
particles. In the course of the transformation virtual pro- 
cesses creating and annihilating quasi-particles are more 
and more suppressed so that the average concentration of 
quasi-particles decreases gradually. Hence towards the end 
of the transformation the neglect of 3- (and more) parti- 
cle terms is well justified. On the other hand, the 3-(and 
more) particle terms are not present in the initial hamil- 
tonian so that their neglect is well justified during 
the first phase of the transformation as well. Only during 
the intermediate phase there is no general control of the 
quality of the approximation. Here one has to focus on the 
particular system under study. To assess the quality of the 
approximation one may either compare to otherwise avail- 
able data on the system (external quality control) or one 
may compute some of the terms neglect in order estimate 
how large they are (internal or self-consistent quality con- 
trol) . In the present work we will use the first method and 
assess the quality of our findings by comparing them to 
the known exact results. Summarizing, the approximation 
is well justified if the system can be considered a dilute 
gas of quasi-particles independent of the actual interaction 
between the quasi-particles. 

Finally, we wish to point out an analogy between our 
approach and a more standard diagrammatic one. By re- 
stricting ourselves to the terms in Eqs. 



dealing with a 2-particle irreducible vertex function which 
is renormalized continuously. The renormalization equa- 
tion (see appendix B) is bilinear in all coefficients. The 
terms contributing to c^-T, which are bilinear in the ver- 
tex function r, result from the commutation and a sin- 
gle contraction of two 2-particle irreducible vertex func- 
tions. In this sense our procedure bears similarities to a 
1-loop renormalization or to the summation of all parquet 
diagrams |37| . The main difference is that our approach 
keeps the problem local in time along the renormalization 
flow because the transformation is unitary. Hence not fre- 
quency dependence enters. We consider this a major ad- 



vantage of the CUT approach for numerical application 
since much less bookkeeping is required since no frequency 
dependence has to be traced. 



3 Numerical Results 

We turn now to the numerical solution of the differential 
equations set up in appendix B. The momentum depen- 
dence of all functions is discretized in an equidistant mesh. 
The positions of the points of this mesh are chosen such 
that no /c-point lies precisely at k = ±fep = ±7r/2. For 
N divisible by 4, this corresponds to antiperiodic bound- 
ary conditions. In this way unnecessary degeneracies are 
avoided. For large system size the influence of the bound- 
ary conditions becomes increasingly unimportant. 

For the system sizes (N w 50) that we will be look- 
ing at there are about 50.000 coupling constants which 
are traced in the set of differential equation. Luckily, the 
differential equations are not very sensitive since they de- 
scribe the convergence to a static fix point. The numerics 
is done by a Runge-Kutta algorithm with adaptive step- 
size control. This algorithm is robust but quite laborious 
because the differential equations have to be evaluated six 
times for each step. Rigorously, the fix point is reached 
at I = oo. In practise, we stop the flow when the relative 
change of E and et falls below 10~ 6 . 



3.1 Illustration of the CUTs 

Fig. [5] illustrates what the continuous unitary transfor- 
mations are doing. The parts of the vertex-function that 
are associated with processes that change the number of 
quasi-particles are transformed away as I goes to oo while 
the other part is renormalized. At the end of the transfor- 
mation there are discontinuities at the borderlines of the 
different parts and some singular kinks if the scattering 
processes occur at the Fermi points or if p — q holds. 

A more quantitative insight is given by Fig. where 
one-dimensional sections of the vertex-function Fk qp are 
shown. 



Fig. 3(a) depicts r 0qq for various values of £. These am- 
plitudes are related to processes where two fermions at 
q and — q are first annihilated and then created again. 
The shape of the function close to \q\ = ir/2 can be fit- 
ted by r 0qq « -0.007q 2 - 0.266 In | \q\/n — 0.5j — 0.125. 
So a logarithmic singularity at ^ can be presumed. 



In Fig. 3(b) the amplitudes of processes are shown 
where two fermions are taken from q and — q and put at 
ir — q and — ir — q. All these processes change the quasi- 
particle number and, hence, have to vanish for £ — > oo. 
The processes near the Fermi wave vector q « ^ are 
decreasing very slowly. 



— For the plots in Fig. 3(c) 



we fix one shift of a fermion 
from below the Fermi level to an energy above it. The 
shift of the other fermion is varied. Here, different 
changes in the number of quasi-particles are possible. 
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Fig. 2. Sections of the vertex function Fk qp at k — f as function of q,p 6 {— tv, tv} at three values of £ (0, 0.3 and oo) for V = 0.5 
and N — 48 sites. The lightest parts are those that change the number of quasi-particles by ±4, the gray parts change them by 
±2, the black parts leave it unchanged. The quasi-particle-number changing amplitudes are transformed to zero. 
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Fig. 3. Sections of the vertex-function Fk qv as in Fk qv : c i+ q c k-q c k- P c k+ P '■■ The pictograms illustrate the scattering processes. 
They show the dispersion (solid line), filled dots for the fermions to be annihilated and open dots for the holes the fermions 
are nut. The arrows denict how the nrocesses change on varvine' the momentum. AOP denotes the change of the numher of 
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1/N 

Fig. 4. Interaction values above which the approximation 
breaks down because of loss of convergence. Filled symbols 
stand for numerically stable runs, open symbols (N > 52) for 
runs where numerical inaccuracies spoil the convergence. The 
dashed line extrapolates the numerically reliable results (see 
Eq.ETJ. 



There are regions where the amplitudes vanish and 
others where they are only renormalized. The process 



at k = 



is apparently different. It corresponds to 



the exchange of the two particles. At present, we can- 
not judge whether this difference will be relevant in 
the thermodynamic limit. It is an interesting question 
whether it retains a finite weight for N — > oo. 



In Fig. 3(d) we show the evolution of a scattering event 
from pure exchange at q = at the right Fermi point 
over forward-scattering for small positive values of q 
to umklapp-scattering at q — it. One realizes that the 
scattering amplitudes at q = and at q — n are partic- 
ularly enhanced by the renormalization. It seems to- 
tally inappropriate to approximate the scattering am- 
plitudes as being constant for small momentum (q « 0) 
or for large momentum (q ss 7r) because the points 
q = and q = ir appear to be decisively different 
from the scattering in their vicinity. This observation 
makes the question whether a continuum description is 
quantitatively applicable an interesting issue for future 
studies. 



is decreasing as N is growing. In addition, numerical prob- 
lems occur at system sizes N > 52 and V close to V C iit.(N) 
due to accumulated inaccuracies. This can be seen in the 
unsystematical behaviour of V CI it.(N) above N = 52 and 
in the unphysical shape of e.g. the dispersion for V close 
to V^ r it.(iV) in Fig.0 On the other hand, the decrease of 
V cr it.(N) is very systematic as long as N < 52 holds. It 
can be approximated very well by a linear dependence 



V ait .(N) w 9.06/7V + 0.75 



(31) 



which leads to an extrapolated value in the thermody- 
namic limit N — > oo of V^ r it.(oo) « 0.75. We conclude 
that it is possible to study the metallic phase of the sys- 
tem and that the conclusions drawn from the finite-size 
calculations are also relevant for the thermodynamic limit. 
A description of the transition to the insultating phase is 
presently not possible. 



The Case N = 4 Some insight on the cause for the loss 
of convergence can be obtained in the simple case of four 
points in momentum space, which is analytically solvable 
with and without approximation. The analytic solution 
is simple because the quasi-particle vacuum \ipo) is con- 
nected by the interaction only to the state where both 
fcrmions are excited \ip e x.}- Thus one has to diagonalize 
a 2 x 2 matrix which can be done easily directly or using 
the CUT. 

Our approximation, however, is dealing with operators 
and not with matrix elements. The state \ipex.) involves 
four quasi-particles - two holes and two excited particles. 
In the approximation 3- and 4-particle terms are neglected 
so that the energy of the excited state is incorrect. Indeed, 
the gap between the ground state and the excited state 
is under-estimated 2 . It even closes for V > y/8 and the 
couplings diverge before reaching I = oo. Thereby, it is 
shown that the approximation may spoil the applicability 
of the approach. 



3.3 Ground State Energy per Site 



3.2 Convergence 

General statements on the convergence of the CUT ap- 
proach chosen are possible 20,2]. But they do not ensure 
that the method works for macroscopically large systems. 
Additionally, the unavoidable use of approximations may 
lead to the break down of convergence. So we investigated 
empirically up to which value of the interaction V^ r it the 
CUT works. This value depends on the system size N or 
the density of the discretization mesh, respectively. If V is 
larger than V CI it.{N) all couplings are diverging at a cer- 
tain I and the solution of the differential equations cannot 
be continued asymptotically to oo. Furthermore, V^ r it.(iV) 



In Fig. 13.31 we plot the shape of the non-linear lowering 
of the ground state energy, i.e. the correlation part of 
the ground state energy, for various system sizes N as 
obtained by CUT. The ground state energy diverges to 
negative values close to V CI it.(N). For V smaller than the 
extrapolated value V cr n.(oo) all systems show nearly the 
same dependence on the interaction V. 

The case N = 48 is compared to the exactly known 
thermodynamic result (271128] in Fig. 13.31 The quantita- 
tive results are very close to each other for quite a large 
region of V. Even for V = 0.5 the relative difference is less 
than 1%. Only as V approaches V crit .(48) = 0.9361 the 
approximate ground state energy is diverging very fast to 
— oo. 
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Fig. 5. CUT results for the correlation part, i.e. beyond linear 
order, of the ground state energy per site for various system 
sizes N. 



Fig. 8. Comparison of the non-linear correlation part of the 
Fermi velocity as obtained for N = 48 by the CUT with the 
exact thermodynamic result, the exact quadratic order V 2 and 
the result obtained from bosonization. 
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Fig. 6. Comparison of the CUT results (N = 48) for the 
correlation part of the ground state energy per site with the 
exact thermodynamic result. 



to 
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k[7t] 

Fig. 7. Renormalized dispersions from the CUT calculation for 
N = 48. The results for V = 0.03, 0.06, 0.09 . . . , 0.9 are shifted 
with respect to each other in order to yield a three-dimensional 
view on the evolution of the dispersion. 



3.4 Dispersion and Renormalized Fermi Velocity 

At the end of the transformation, i.e. at I = oo, when 
the effective hamiltonian has reached block diagonality 
only scattering processes are left which leave the num- 
ber of quasi-particles unchanged. Then Sk is the renor- 
malized 1-particle dispersion. This means that in the ef- 
fective model after the transformation it is possible to 
add a single quasi-particle (hole or particle) to the ground 
state such that the resulting state is an exact eigen state 
because there is no other quasi-particle to interact with. 

Note that this statement is not in contradiction with 
the widely known fact that the single-particle propaga- 
tor G(k,u) of a Luttinger liquid does not display quasi- 
particle peaks p^ll38U89lfT7] because the single-particle prop- 
agator G(k,oj) refers to adding or taking out a fermion 
before any transformation. It is an interesting issue, yet 
beyond the scope of the present work, to apply the CUT 
to the creation and annihilation operators in order to 
recover the usual Luttinger liquid result in the framework 
of the CUT renormalization. 

Fig. [7| shows the evolution of the dispersion for N = 
48 on increasing V. The dispersion behaves like a cosine- 
function with a renormalized Fermi velocity as expected 
from the exact result. For V close to V^ r it. (48) (at about 
V = 0.7) there are kinks emerging. We reckon that these 
kinks represent spurious features induced by accumulated 
numerical inaccuracies for the same reasons for which the 
convergence is hampered for large system sizes. 

Renormalized Fermi Velocity By fitting the function 
—Up cos(fc) to the calculated dispersions we obtain results 
for the renormalized Fermi velocity v F as function of the 
interaction V. These values are dominated by the linear 
Fock term 1 + — or its finite-size equivalent in Q14|). re- 
spectively. In order to yield a better resolution of the influ- 
ence of correlations we substract the constant and linear 



2 But the deviation occurs only in fourth order in V. 
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Fig. 9. Momentum distribution as obtained by CUT for V = 
0.6 and 48 sites. Lines are interpolations to the three points 
closest to the Fermi wave vector assuming Luttinger l|33fl or 
Fermi behaviour (I34L respectively. 
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Fig. 10. Luttinger exponent a(V) as obtained from least- 
square fits of the CUT momentum distributions. 



3.5 Momentum Distribution 

Last but not least we analyze the momentum distribu- 
tion n(k) in order to show that Luttinger liquid behaviour 
is retrieved. The standard approach to do so would be 
to transform the operator c\c k according to Eq. (|3J) ^3 
140) . Here, however, we use a simpler approach suitable 
for static expectation values and correlations. The value 
Eoo — E{1 — oo) is the ground state energy within our 
approximation because the block diagonal effective hamil- 
tonian H(£ — oo) does not contain processes exciting the 
vacuum any more. First order perturbation theory shows 
straightforwardly that functional derivation of yields 
the expectation value n(k) 



rf£oo = /dH\ 
de k \ de k / 



:) =: n(k) 



(32) 



In the numerical treatment the functional derivation can 
be easily realized by approximating the ratio of infinites- 
imal differences by the ratio of small finite differences. 
Hence no serious extension of the algorithm is needed to 
compute the momentum distribution. For relative varia- 
tions of Ek around 1CU 3 the ground state energy is lin- 
ear in these variations within an accuracy of about 10~ 5 . 
So n(k) can be determined to this accuracy. In Fig. ED 
the momentum distribution (symbols) for N = 48 and 
V = 0.6 obtained in this way is depicted. 

Due to the discretization n(k) can be computed only 
for a finite set of points. At first sight, a jump seems to 
dominate at the Fermi wave vector tt/2. But a discretized 
power law distribution displays also a jump - in particular 
if the exponent is small. In order to understand the na- 
ture of the distribution a quantitative analysis is required. 
Hence we fit the distribution obtained to two functions, 
one being appropriate for describing a Luttinger liquid [321 

El 



\n{k) 



C\(Ak) a +C 2 Ak 



(33) 



the other being appropriate for describing a Fermi liquid 

EH 



terms. The remainder v F _ rm is plotted in Fig. [5] It is com- 
pared to the exact result the exact quadratic term V 2 
and the result obtained from bosonization (cf. appendix 
C). The shape of the curves is basically the same. The 
CUT result is too small (in modulus) compared to the ex- 
act results. For small V this is mainly due to a finite-size 
effect as comes out from an extrapolation N — > oo. Again, 
the CUT approach is less reliable close to the critical in- 
teraction value V CI it(N). 

Note that the bosonization results fits less well to the 
exact result than does the CUT result. This is due to 
the fact that in bosonization only the processes infinitely 
close to the Fermi points are considered. The deviation 
between the dash-dotted second order curve from the ex- 
act results reveals that higher order terms are important 
as well. They are partly captured by the CUT procedure. 



\n(k) 



C 1 \n{Ak)Ak + C 2 Ak (34) 



with Ak := \k ~ fcp|. For both possibilities three free pa- 
rameters (a, Ci, C 2 or Zk F , C\, C 2 , respectively) are de- 
termined. In Fig.Elthe parameters are fixed to interpolate 
the three points closest to the Fermi wave vector. Another 
way is to perform a least-square fit; the resulting curves 
are shown in Fig. 3 in Ref. JSj- In both analyses, the qual- 
itative result is the same. The Luttinger fit describes our 
data much better than the Fermi fit. We conclude that 
our data describes rather a power law behaviour than a 
jump. We cannot exclude, however, a behaviour compris- 
ing a power law behaviour and a jump. But there is no 
reason to believe that such a behaviour should occur. 

Due to the small exponents a occurring and the re- 
stricted system sizes the power law cannot be distinguished 
reliably from a logarithmic behaviour or from a function 
of some logarithm in |fe — kp\. Note, however, that the 
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simple logarithm found in Ref. 10 is not likely to occur 
since we do not transform the observable in leading order 
only. Infinite orders of the interaction V contribute to the 
ground state energy and hence to the derivative 

In order to push the analysis one step further we com- 
pare the exponent a resulting from the fits to the exact 
one. The points in the vicinity of the Fermi wave vector 
are more influenced by finite size effects 02] ■ Eventually, 
we choose the exponents coming from least-square fits. In 
Fig. 1101 thev are compared to the exact values, to the sec- 
ond order result and to the values coming from the direct 
application of bosonization (cf. appendix C). Clearly, the 
CUT results agree very well with exact data for not too 
large values of V < 0.6. From the comparison to the exact 
second order term V 2 one sees that the CUT data de- 
scribes the full exact result better than the second order 
term alone. We infer that the CUT data reproduces the 
third order term also. This does not come as a surprise 
since we showed above that the ground state energy is 
exact including the V 3 term. Thus any quantity derived 
from it will also be exact in the same order. So the mo- 
mentum distribution n(k) and hence the exponent a have 
to be exact up to and including V 3 . 

The comparison to the bosonization result in Fig. I1UI 
is also instructive. By construction the bosonization result 
captures only the physics in an infinitesimal vicinity of the 
Fermi points. There no umklapp scattering is possible as 
can be seen from Eq. (|13(l or from the analyses in Rcfs. 

(cf. discussion after Eq. JTEJl). For this reason, the 
bosonization fails to detect any precursor of the incipi- 
ent phase transition to the CDW occurring at V = 1. 
The dependence of the a obtained from bosonization is 
hence much smoother than the exact result. The CUT re- 
sult includes scattering at all momenta. So precursors of 
the phase transition can be captured. Unfortunately, the 
breakdown of the approximation as used here makes fur- 
ther statements on the description of the phase transition 
impossible. 



4 Discussion 
4.1 Conclusions 

We presented general arguments in favour for a renor- 
malization treatment of interacting fermionic systems by 
means of a continuous unitary transformation (CUT). The 
precise choice of the unitary transformation, i.e. the choice 
of its infinitesimal generator, was motivated and the gen- 
eral properties of the transformation were elucidated. The 
method was illustrated for a model of one-dimensional, re- 
pulsively interacting fermions without spin at half-filling. 
The technical considerations for the explicit calculation 
were given. An important point was a consistent, non- 
redundant notation l|18|) which made the crucial cancella- 
tion between the Cooper pair and the zero sound channel 
manifest. Results were obtained for the correlation part of 
the ground state energy, for the 1-particle dispersion, for 
the 2-particle vertex function and for the static momen- 



tum distribution. The findings were compared to exact 
results as far as possible. The agreement was very good. 

The CUT employed uses standard quasi-particles as el- 
ementary excitations. So it represents an explicit construc- 
tion of Landau's mapping of the non-interacting quasi- 
particles to the elementary excitations of the interacting 
system jT^j. Note that the existence of such a smooth 
connection is not surprising since already the bosoniza- 
tion identity for fermionic field operators fMT7\ repre- 
sents such a smooth link. In the continuum limit, a spin- 
less model with or without interaction can be mapped to 
a single mode boson model with linear dispersion. Hence, 
it is possible to link the interacting model via the bosonic 
model to the non-interacting one (see also the remark on 
Kchrein's results |48II44| below) . We take the success of our 
approach as corroborating (numeric) evidence that Lan- 
dau's mapping exists in one dimension. In particular, the 
numerical findings of the momentum distribution n(k) in- 
dicate that important features of Luttinger liquids could 
be retrieved. Signatures of a Luttinger-type power law be- 
haviour at the Fermi points were found, though hampered 
by the accessible restricted system sizes. 

We are convinced that CUTs represent a powerful renor- 
malization scheme for low-dimensional systems. Since no 
states are eliminated the effective hamiltonian obtained at 
the end of the CUT allows to compute spatial (shown here) 
and temporal correlations (for an example, see Ref. |14j ) 
at small and at large wave vectors or excitation energies, 
respectively. So, in principle, no information is lost in the 
course of the renormalization, in contrast to, for instance, 
Wilson's renormalization Qj. Of course, approximations 
which are necessary in practical calculation will introduce 
some uncertainties. Recent developments in renormaliza- 
tion approaches by integrating out degrees of freedom al- 
low also to compute high energy features if the observ- 
ables are equally subject to the flow, see for instance the 
appendix in Ref. A major advantage of the CUT ap- 
proach is that no frequency dependence needs to be kept. 
This represents an important facilitation for the actual 
numerical realization of the renormalization. 



4.2 Connections to other Work 

Life-Time of Excitations If a complete or partial diago- 
nalization is obtained by a continuous unitary transforma- 
tion the eigen values are by construction real. So the exci- 
tations are well-defined in energy and they do not display 
a finite life-time. This statement appears almost trivial if 
one bears in mind only the mathematical linear algebra. 
From the physics points of view, however, one might be 
surprised since one is used to that excitations, for instance 
quasi-particles, have a finite life-time in many-body sys- 
tems. This comes about because no true eigen state are 
considered when an excitation of finite life-time is studied. 

For instance, adding a fermion to the ground state of 
an interacting fermion system by application of a sim- 
ple creation operator generically does not yield an eigen 
state but a sophisticated superposition of true eigen states. 
This looks as if the "true" fermion added decayed because 
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its spectral function displays a peak of some finite width. 
Technically, one still uses the normal single-particle ba- 
sis but the self-energy E(k,oS) acquires a finite imaginary 
part which represents the fact that the single fermion is 
coupled to states with one or more particle-hole pairs. One 
may view + S{k, uj) as an imaginary eigen value of the 
single particle state. But the single particle state is not an 
eigen state of the underlying hamiltonian. 

The description using a unitary transformation de- 
signed to diagonalize the original hamiltonian is different 
from the standard approach in physics sketched above. 
One tries to find real eigen vectors and eigen values. Broad 
spectral functions do not occur because of imaginary eigen 
values but because of the superposition of eigen vectors. 
Kehrein and Mielke coined the expression that "the ob- 
servables decay" and described the phenomenon in the 
context of dissipation |4fi| . 

So it is not surprising that, after the CUT is applied 
to one-dimensional fcrmions, we find quasi-particles with 
infinite life-time. These are quasi-particles after the trans- 
formation. The fermion before the transformation will de- 
cay into states with additional excited particle-hole states. 
So there is no contradiction. In this context it is interesting 
to note that there is a formulation of diagrammatic per- 
turbation theory which uses also infinite life-time excita- 
tions |47U48| . In this formulation the fermionic one-particle 
states acquire a renormalized eigen energy due to the in- 
teraction. The change of the eigen energy e — ► e' leads 
to a change in the occupation number (1 — exp(— [3{e — 
A 1 ))) -1 -* (1 - exp(-/3(e' - m))) -1 ) 47 . So this approach 
corroborates that a description of the physics in terms of 
non-decaying quasi-particles is possible. The question how 
dynamic correlations can be described is not discussed in 

Refs. mm. 



CUTs as Numerical Approach While the present work 
was being finished White ^5] proposed a numerical scheme 
which is very similar in spirit to what we did here. Besides 
discrete transformations which work less efficiently, a con- 
tinuous unitary transformation is used with an infinitesi- 
mal generator with matrix elements 



1 



Hij(£) 



(35) 



where a 1-particle basis is used, i.e. a basis in which the 1- 
particle part of the hamiltonian is diagonal. The 1-particle 
eigen energies are given by Ei. The approach is applied 
to a small molecule, namely H2O, and the ground state 
energy is computed very reliably by rotating the original 
ground state to a Fermi sea, i.e. the quasi-particle vacuum. 
That means that states which partially occupied n > 0.5 
are mapped to filled states and states which are partially 
empty are mapped to empty states. This is what we did 
in our present work as well. So Ref. |15| provides an inde- 
pendent investigation of the power of continuous unitary 
transformations for a different fermionic system. In addi- 
tion, it is investigated in Ref. to eliminate a number of 
states which lie far off the Fermi level without diagonaliz- 



ing the problem completely by the continuous transforma- 
tion. The remaining effective problem, which is simplified 
considerably due to the reduction of the Hilbert space, is 
then solved by standard diagonalization algorithms, e.g. 
DMRG. Also this approach proved to be very powerful. 



Choice of CUT In Ref. JU|, Wegner investigated a one- 
dimensional n-orbital model in the continuum limit by a 
continuous unitary transformation. The main difference 
in the unitary transformation is the use of a generator 
different from the one in Eq. ©, namely 



(36) 



where H-q is the part of the Hamilton one wants to keep. 
The approach succeeded when Hd comprised all terms 
that do not change the number of quasi-particles. Hence 
this renormalizing scheme is very similar to the one used in 
the present work. It would be an interesting issue to com- 
pare both approaches quantitatively in a simple model. 
There are arguments in favour for both of the two ap- 
proaches. 

First, the choice © has the advantage that the kind 
of terms that are generated is restricted: the block band 
structure is preserved. Second, off-diagonal terms chang- 
ing the number of quasi-particles arc eliminated even if 
they are not accompanied by a change of the 1-particle 
energies, i.e. certain degeneracies are lifted. Third, a sup- 
pression of off-diagonal parts starts already linearly in the 
differences of the diagonal parts. A weakness of the choice 
© is given when there are scattering processes which in- 
crement the number of quasi-particles but decrease the 
1-particle energies. In this case, the corresponding ampli- 
tudes are first enhanced before the decrease to the end of 
the transformation. Due to the necessary truncations it 
may be difficult to control the quality of the approxima- 
tion during the stage of enhancement. 

On the other side, the choice l|36(l is firstly very ro- 
bust since off-diagonal terms are always suppressed due 
to the fact that the energy difference occurs squared |l(Jj . 
Second, one does not need to know explicitly the eigen 
basis of Hd. The squares are generated automatically by 
the double commutator when Eq. I|36|l is combined with 
Eq. (J2J). Yet these two commutators must also to be com- 
puted which might be tedious. Another weakness arises 
when large-scale degeneracies spoil the method by stop- 
ping the renormalization prematurely. For instance, the 
vanishing of 77 implies the stop of the flow but guarantees 
only that there is a common basis set of H-q and H , not 
that H is diagonal. So the conservation of the number of 
quasi-particles is not ensured by the choice (|3l)|l . 

Summarizing the comparison of the choices H6I36J) we 
reckon that JHJ works better if the number of excitations 
correlates well with the energy. If this is not the case, the 
robustness of H3HJ1 may be preferable. Note that there are 
still completely different generators conceivable For 
example, one can use in a basis of 1-particle states 



rjij (£) = if gi = q 3 
f]ij(£) = sign(£' J - Ej)Hij 



otherwise 



(37) 
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where in case of degeneracy Ei — Ej the energy of the 
state with more excitations is assumed to be mfmitesi- 
mally higher. Such an approach captures advantages of 
Ipjjl while avoiding its disadvantage. Further investigations 
of these issues are certainly called for. 



Fermionic Excitations The fact that we treat a system 
of interacting one-dimensional fermions in its Luttinger 
liquid phase without using explicitly collective bosonic 
modes might be surprising. Yet it is not uncommon that 
an interacting system is suitably described (after certain 
transformations) by free or nearly free fermions. As an 
example we quote the work by Kehrein who succeeded 
to map a sine-Gordon model by a sequence of continu- 
ous transformations onto a model of free fermions |43II44| . 
Indeed, Kehrein's mapping accomplishes this aim for a 
broader range of parameters than previous renormaliza- 
tion treatments. The neglect of interactions between the 
excitations viewed as elementary is a fairly severe approxi- 
mation for 1 < An where the breathers, i.e. bound states, 
are known to occur |50| . Their description requires the in- 
clusion of the interaction between elementary excitations. 



Landau's Fermi Liquid Finally, we wish to comment on 
the use of a Landau's Fermi liquid description in terms 
of quasi-particles for one-dimensional systems |13| . The 
possibility and the power of such a description has been 
noted previously by Carmelo et al. [HI] in the framework of 
the Bethe ansatz solution of the one-dimensional Hubbard 
model. Carmelo and coworkers use the spinons and holons 
as they arise in the Bethe ansatz solution as pseudo-particles 
Then an approximate treatment for the spatial and tem- 
poral |52| correlations is built by describing the excitations 
as small deviations from the ground state distributions of 
these pseudo-particles. In this sense, the concept of Lan- 
dau's Fermi liquid is generalized to one dimension. The 
author emphasize, however, that the pseudo-particles can- 
not be smoothly linked to the gwasi-particles of the non- 
interacting solution. In this point, a clear difference to our 
finding here and in Ref. ^2] occurs. We argue on the basis 
of our numerical results that a smooth mapping between 
the interacting and the non-interacting excitations exists 
even in one dimension. The existence of such a mapping 
as long as the system remains massless can already be 
deduced from bosonization. Since the interacting and the 
non-interacting model can be mapped to a model of free 
linear dispersion bosons (discarding Umklapp scattering) 
they can also be mapped to each other. By means of the 
CUT we constructed such a mapping explicitly. 

The generalized Landau liquid in Refs. {51} relies on 
the Bethe ansatz solution of the one-dimensional Hub- 
bard model. Hence it may be that the integrability is a 
prerequisite for the generalized Landau liquid. In our cal- 
culation in contrast, the integrability of the model studied 
does not play a role other than providing a rigorous bench- 
mark. But so far, we have not considered the spinful case. 
Its investigation by CUTs is certainly called for. 



4.3 Outlook 

A comprehensive summary is given in Sect. 4.1. Here we 
point out in which directions further work is required. In 
view of the numerical nature of the present work, an ana- 
lytical treatment would be helpful. There is still a certain 
gap between the analytical result of a logarithmic diver- 
gence for the momentum distribution obtained in Ref. |10| 
and the numerical results we found. For an analytical 
treatment the models to be considered have to be sim- 
plified further. Spin, however, should be included in order 
to enlarge decisively the class of systems which can be 
described. 

Another very interesting issue is the computation of 
dynamical quantities like the local spectral function A(u) 
or the momentum resolved spectral function A(k,uj). In 
the framework of standard renormalization such investiga- 
tions are presently carried out |53|. The dynamical quan- 
tities are of interest to see theoretically to which extent 
and to which accuracy they can be computed at all en- 
ergies and momenta. For the explanation of experimental 
data the spectral functions are of utmost importance. It 
is this objective which requires in particular to go beyond 
the asymptotic regime of very small energies and momenta 
54 . For gapful spin systems dynamical quantities have al- 
ready been computed successfully [T2|- The- results agree 
very well with experiments and render deeper insight in 
the underlying physics [SSIES] • 

In order to go beyond one-dimensional systems, mod- 
ified generators have to be investigated, see e.g. |15]. The 
pros and cons of the choices used presently were briefly 
discussed in the preceding section. The issue of the opti- 
mum generator represents a longer-lasting question since 
the answer depends certainly on the model to be studied. 
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A Normal-Ordering 

Normal-ordering |37] denoted by colons : A : is a standard 
procedure which is explained at length in the text books 

. So we recall here only the gist of it which is technically 
relevant for our calculation. For further details we refer the 
reader to the concise script by Wegner |59| . 

Considering a product X n of n fermionic operators it 
is a priori not clear on how many particles this operator 
really acts. It looks as if it acted on n fermions. But a 
part of this action may be redundant in the sense that 
it can be expressed also by an operator Y m with m < n. 
This is indeed the generic situation. A point of reference 
is needed in order to be able to define how many particles 
are involved in a certain process. Given the ground state 
of a 1-particle hamiltonian 3 the normal-ordering ensures 
that the normal-ordered n operator does not contain parts 
which can be viewed as action of an operator with less 
fermionic factors. This is the physical content of 

(: P n :: Q m :) = for n ? m (38) 

as derived for normal-ordered terms |59| . 

Technically, a usual product of m fermionic operators 
ak is expressed in terms of normal-ordered terms as 

«feiafc2 • ■ ■ dfem = 

: exp Gu gaI iS aalolt a kl a k2 ...a km :, (39) 

\M 1 1 I 

3 The formalism works identically for finite temperatures 
with respect to the statistical operator of a 1-particle hamilto- 
nian. 
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where G k i is the contraction (a^a;). The superscripts 'left' B The Self-Similar CUTs 

and 'right' indicate that in the double derivatives only 

pairs are taken where a^ lt is a factor to the left of af ght . To In the self-similar approbation we compute 
obtain the correct signs the 'left' derivation must be taken 
before the 'right' excitation. Eq. (|39|l stands for the known 
procedure that a product is normal ordered by writing 
down the sum of terms with all possible numbers and sorts 
of contractions. The inverse relation is given simply by 



+ 2 



+ 2 



: o. k \a k 2 ■ ■ ■ dkm '■— 

eXP ~ H Gkl da'^da^ afclflfc2 ■ ' ■ akm ■ ( 40 ) 

V hi 1 1 J 

Combining Eqs. (|39|l and l|4U|l leads to the useful expres- 
sion for products of normal-ordered terms : A(a) : and 
: B(a) : 



T ( %) +T^\i) +^T 2 ( i 2) W,^sign(z)T 2 ( f ) (f) 

-2 -2 

(43) 

(3) 

neglecting the arising terms. The change in the num- 
ber of quasi-particles for the expression 

is given by S kqp as defined by 



r, 



A(a) :: B(a) :=: exp 




A(a)B(b) 



b—a 

( 41 ) 

where the superscripts 'left' and 'right' are no longer needed 
due to the sequence of factors in the product. In practice, 
Eq. I|4I|) means that for the normal-ordering of a product 
of already normal-ordered factors not all contractions need 
to be considered. Only those contractions matter where 
the two fermionic operator do not come from the same 
factor. This is easy to understand since the contractions 
between fermionic operators from the same factor are al- 
ready accounted for by the normal-ordering of each factor 
separately. 

In order to determine the differential equations result- 
ing from Eq. (|2J) commutators of normal-ordered terms 
must be computed. To do this without passing by non- 
normal-ordered expressions an extension of Eq. (|41|l to 
commutators is particularly useful. We derived and checked 
the identity 



A{a) :,:B(a) :] = 

:eX P {^ Gkl ( da^db^ + 

[A(a),B(b)} 



ht p. loft 



db" sM da 



b—a 



(42) 



where the commutator is computed using the anticom- 
mutators {a k ,h} '■— {o-k,o,i}. Eq. (|42|l means that one 
can first compute the commutator as usual, but remem- 
bering whether the fermionic operator comes from A or 
from B. Then normal-ordering is achieved by writing down 
the terms with all possible contractions between pairs of 
fermionic operators where one comes from A and the other 
from B. In this way, the computation of the actual gen- 
eral flow equation J2J becomes a task which is not too 
demanding. 



kqp '■ c k+q C k-q C k-p C k+p 



Skqp = sign(n fc+p + - n k - q - n k+q ) , (44) 

where we use n k for the momentum distribution of the 
unperturbed Fermi sea since this is the quasi-particle vac- 
uum to which we are mapping the ground state. (The ac- 
tual momentum distribution is denoted n{k)). With these 
definitions we calculate the commutator using normal- 
ordering as explained in appendix A. Then we compare 
the coefficients of the various terms (O-particle, 1-particle 
and 2-particle terms) and determine in this way the set of 
differential equations. 



B.l Ground state energy per site 

The differential equation of the ground state energy per 

(2) (2) 

site depends only on the commutator [T± 4 , T^ 4 ] where all 
fermionic operators are contracted. There are four differ- 
ent ways to combine the operators for the complete con- 
tractions. But all of them lead to the same expression if 
one uses the symmetries Ea. (|15f) - (I18|I . Finally one obtains 



d E 

diN = Tf3 ( 1 - 2n k+p)nk+ q n k - q S k ^ p rl qp .(A5) 



N 3 Z/ 

fce[0,7r) 



B.2 Dispersion 

For the 1-particle term, the dispersion, one has to take 
all combinations of three contractions into account that 
occur in [T^ 2 ) 4 , T^} with i e {0, ±2, ±4}. One obtains 16 
different parts that turn out to be identical. In order to 
avoid double-counting one of the free momenta must be 
restricted to [0, 7r) or the sum must be divided by two. So 
one obtains finally 



dl £k 



Jp Y ~ 2ri k+q ^ p )n k - 2 q 

qp&[— 71", tt) 



Sk- 



qqp 1 k-qqp 



Tln+p—qTln—p—q) X 

(46) 



B.3 Vertex Function 



For the vertex function B kqv there are two commutators to 



calculate: [if 2 } ,Tf } ] (i, j e {0,±2,±4}) with all possible 
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combinations of two contractions and [Tq , T^] with all 
possibilities of a single contraction. To keep the notation 
short, we define 

^kqp \kqp ■— (Skqp — SKQp)rk qp r K Qp . (47) 

The differential equation for the flow of F kqv then reads 

~^J^kqp — (&k-\-p H~ ^k—p &k — q ^-k+q^^kqp-^kqp ~t~ 



Qe[-7r,7r) 



n - ,j | n r« — 
2 2 


W K — « — 

2 




-Pi 


" 7 

2 


2 


' . 


-g 


r. ■ p 

2 




— q + Q k — q- 


Q fc + q- 


Q 


+p| 


fc-p+Q 


fc + p- 


Q 


+q 




Q 


2 2 


2 




2 


2 




2 




fc + q + Q fc + q- 


Q fc-q- 


Q 


+P| 


fc-p+Q 


fc+p- 


Q 


-q 




Q 


2 2 


2 




2 


2 




2 




k — q + Q k — q- 
r> 


Q fc + q- 
9 


Q 


-/-'I 


k+P+Q 

9 


9 


Q 


+9 


9 


Q 



} + 2(1 - 2n Q )$ kq{Q - kmQ _ k)p . (48) 

Note that the appearance of four <£-terms is due to the fact 
that we denote the scattering processes in the hamiltonian 
(TjQ) in a notation symmetric way H18|) . A naive compari- 
son of coefficients would lead only to one of the four terms. 
The other three come into play if one requires that 7§Tj. gj , 
fulfills l|18fl . As explained in the main text, the notation 
obeying (|18fl ensures that a maximum number of cancel- 
lations are dealt with explicitly. This is advantageous on 
the numerical as well as on the conceptual level. 



C Bosonization of the Model 

We employ a constructive bosonization by linearizing the 
dispersion on both branches (r = +1 «-» right branch, 
r = — 1 «-> left branch) 

£k,r=r(l + ^)(k-rk F ) . (49) 

Furthermore, the bare interaction vertex (|13fl is evaluated 
with all momenta being taken at the Fermi points ±k-p. 
The corresponding scattering strengths take the value ± V. 
Following Ref. |17) one can then determine the renormal- 
ized Fermi velocity as 

4(F) = (l + ^l-(^) 2 (50) 

and the exponent a occurring in the momentum distribu- 
tion from Eq. l|2*2|l and 

, « 1 / 2tt , s 

= 2 V^TW • (51) 



